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ABSTRACT 

Laminar boundary layer separation phenomena in oscil- 
lating flow was experimentally investigated. Multiple hot- 
wire anemometers were used to obtain instantaneous boundary 
layer velocity profiles on a model in an oscillating 
freestream. Certain instantaneous profile behavior was 
found to be uniquely related to wake formation, while non- 
wake-inducing transient flow reversals were found to occur 
throughout the adverse pressure gradient regime. Based on 
these observations, a practical definition of boundary 
layer separation in oscillating flows was formulated. 
Separation data then obtained revealed that the imposition 
of oscillations caused the separation point to move upstream 
of its steady flow location. From the data acquired, a 
parametric correlation of separation point location was 
developed. 
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I. 



UNSTEADY LAMINAR BOUNDARY LAYER FLOWS 



A. INTRODUCTION 

Knowledge of the character and behavior of the boundary 
layer in unsteady flow is of prime importance to the full 
understanding of many practical problems in fluid mechanics. 
Historically, rotating stall in turbo-machinery, aero- 
dynamic flutter, and the flow around accelerating bodies have 
been the most prominent of the problems in this class. 

Recent increases in use and rapid development of rotary wing 
aircraft, with their inherently complex aerodynamics, have 
re-emphasized the need for such knowledge. 

B. PREVIOUS ANALYTICAL STUDIES 

Most of the significant contributions to the resolution 
of these problems are analytical investigations into the 
response of boundary layers to large scale freestream 
unsteadiness resulting from accelerating or periodically 
fluctuating flows. The earliest analysis of such flows 
is the classic work of Stokes [1] who treated an infinite 
lamina oscillating in its own plane in a fluid at rest. 
Rayleigh [2] considered the second order effects in this 
same problem. Schlichting [3] expanded Rayleigh's analysis 
by employing boundary layer assumptions. These early 
analyses treated what is, in essence, the simplest special 
case of the general problem, unsteady viscous flow in the 
absence of a mean flow or pressure gradients. 
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The first significant analysis of unsteady boundary 
layer flows that included the presence of a mean flow was 
that of Lighthill [4] who treated the case of small, low 
frequency sinusoidal oscillations superimposed on the mean 
flows about a cylindrical body, and a semi-infinite flat 
plate. His analysis was based on a "small perturbation" 
treatment of the boundary layer equations for sinusoidal 
flow, retaining only first order perturbation terms. For 
the lower ranges of frequency, a solution was obtained using 
a Karman-Pohlhausen method. For higher frequencies, the 
relatively small convective terms were ignored, resulting 
in equations similar to those of Stokes' "zero mean flow" 
solution. A nxamber of fundamentally important results were 
determined; at low frequencies, the boundary layer flow 
is essentially quasi-steady, having the same character- 
istics for any instantaneous magnitude of freestream 
velocity as it would have in a steady flow of the same 
velocity, and is independent of the dynamic history of the 
freestream flow. As the frequency of oscillation is 
increased, a limiting frequency is reached, beyond which 
the boundary layer responds to the oscillations as it 
would in the absence of a mean flow. This limiting solution 
is commonly known as the "shear wave" solution. Lighthill 
found that the maximum value of skin friction at any point 
leads the maximum in freestream velocity, presumably due to 
a faster response of the lower velocity laminae in the 
boundary layer to changes in the time-dependent pressure 
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gradient than the response of the freestream to these 
changes. Lighthill’s work has been treated as a point of 
departure for many subsequent investigations. Lin [5] has 
analyzed the high frequency case for large amplitude 
oscillations and confirmed the existence of the shear 
wave solution. He also concluded that the governing 
equations become essentially linear at high frequency, 
allowing separation of the effects of the oscillatory and 
steady components of the flow. Carrier and Di Prima [6], 
using an analysis based on a system of equations derived 
by a modified Oseen linearization of the equations of 
motion, confirmed both the shear wave solution and the 
phase advance of wall shear found by Ligh thill. Nickerson 
[7] has expanded the scope of Lighthill's flat plate 
analysis by initially retaining perturbation terms beyond 
first order. Quite recently, Barriol and Lucius [8] have 
used numerical methods to obtain asymptotic solutions to 
the boundary layer equations for oscillating flow, remote 
from the leading edge of a semi-infinite flat plate, in 
the absence of a pressure gradient. Their solution is in 
close agreement with the results of Lighthill and Nickerson, 
but has a much broader range of applicability. 

The response of the laminar boundary layer to small 
freestream oscillations in the vicinity of a stagnation 
point has been investigated by both Rott [9] and Glauert 
[10]. The exact solution for a stagnation flow is 
demonstrated by Glauert to be obtainable from a single 
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ordinary differential equation. Glauert's solutions for 
high and low frequency limiting cases, as well as skin 
friction phase angle lead, agree with Lighthill's 
calculations . 

Compressibility effects in the laminar boundary layer 
for nearly quasi-steady flow were considered by Ostrach 
[11] , who treated first order deviations from quasi-steady 
characteristics. His work again confirmed the skin 
friction phase angle lead found by others. Moore [12] 
improved the analysis by including higher order effects 
for both isothermal and adiabatic flat plates. He obtained 
time-average velocity profiles, in close agreement with 
Blasius profiles, for the quasi-steady regime. Very 
recently King [13] has obtained similarity solutions for 
low speed compressible oscillating flow and for the 
special case of an oscillating wedge in a hypersonic flow. 

Koob and Abbott [14] have recently developed a unique 
integral method for the analysis of the time dependent 
laminar boundary layer in incompressible flow. The method 
is demonstrated for the case of an accelerating flat plate, 
and shows promise of being extendable to more complex 
unsteady flow problems 

C. PREVIOUS EXPERIMENTAL STUDIES 

In view of the rather limited scope of analytic studies 
of unsteady boundary layer phenomena, it is somewhat sur- 
prising that more experimental investigations have not been 
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attempted. The relative lack of experimental data is 
undoubtedly due to the problems encountered in establishing 
the proper flow environment and accurately recording the 
resulting phenomena. 

Richardson and Tyler [15] conducted one of the earliest 
experimental studies in confirming the analyses of Grace 
[16] and Sexl [17], who treated the problem of oscillating 
flow in a tube in the absence of a mean flow. Nickerson 
[7], in conjunction with his analytical work, made hot-wire 
measurements of the laminar boundary layer on a flat plate 
oscillating in a Blasius type mean flow and obtained time- 
average velocity profiles and heat transfer data for the 
quasi-steady (low frequency) regime. His work confirmed 
in part, both his and Lighthill's analyses. Mechanical 
difficulties precluded a more substantial experimental 
program. This work was subsequently extended by Hill [18] 
who imposed oscillations directly on both Blasius and 
Howarth type mean flows. Hill employed a tubular wind 
tunnel equipped with an oscillating gate valve, capable 
of producing frequencies of up to four c.p.s. In the case 
of the Howarth mean flow, the pressure gradient was con- 
trolled by a porous tube along the axis of the test section. 
Hill's work with the Blasius mean flow confirmed both 
Nickerson and Lighthill's analytical results. 

Hori [19] has investigated the laminar boundary layer 
on a circular cylinder oscillating at low frequency in 
steady flows of extremely low velocity. His results 
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confirm Lighthill's quasi-steady regime solutions for 
oscillating cylinders. 

D. TRANSITION AND SEPARATION 

Relatively little work has been reported on two very 
important aspects of unsteady boundary layer flow, 
transition and separation. Many investigations have been 
conducted on problems involving these phenomena, such as 
"retreating blade" stall, rotating stall, flutter, and 
airfoil performance in the presence of turbulence, but 
they have in most cases been restricted to gross effects, 
and rarely treat boundary layer phenomena. 

Miller [20] investigated transition phenomena on a 
flat plate in a Blasius mean flow on which oscillations 
were imposed by using a rotating shutter valve. The 
design of the shutter valve allowed Miller to attain 
oscillation frequencies from 4 to 125 c.p.s. and ampli- 
tudes from 8 to 67 per cent of the freestream velocity. 

This broad range of operating conditions makes this work 
of particular interest and significance, as prior experi- 
mental efforts had been mechanically restricted to 
frequencies of less than 20 c.p.s. and amplitudes no greater 
than 10 per cent of the mean freestream velocity. Miller 
confirms the analytical predictions of Liepmann [21] and 
Greenspan and Benney [22]. In addition, he has determined 
the Transition Reynolds Number, Turbulent Reynolds Number 
and turbulent intermittency factor for oscillating Blasius 
type flow. 
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Even less information exists on separation phenomena 
in unsteady flow than on transition. Analytical approaches 
of proven merit for oscillating Blasius flows have in 
general proved to be inadequate for the treatment of 
oscillating flows in the presence of adverse mean pressure 
gradients. Lighthill's [4] analysis of the laminar 
boundary layer on a cylindrical body in oscillating flow 
loses its applicability soon after the pressure gradient 
on the body becomes positive. Similarly, Hori ' s [19] 
analytical predictions for a cylinder in oscillating flow 
diverge markedly from his experimental results as soon as 
adverse pressure gradient effects are encountered. Hill 
[18] has analytically treated the case of oscillations 
superimposed on a Howarth mean flow and experimentally 
confirmed his solution for one relatively small adverse 
pressure gradient. He does not, however, consider his 
analysis applicable in the neighborhood of separation. 

Both Hill and Hori, although conducting experimental 
investigations in adverse pressure gradient regimes, 
present no data relevant to laminar boundary layer separ- 
ation. It should be noted, however, that these experiments 
were conducted in conjunction with analytical investigations 
and were therefore restricted to the limits of applicability 
of the analyses. All of the investigations mentioned above 
•have been primarily concerned with determining changes 
occurring in the maximum amplitude and phase relation of 
the oscillating velocity component in the boundary layer. 
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It would appear that knowledge of separation phenomena can 
most readily be acquired by observing the changes occurring 
in the complete velocity profile during a full cycle of 
oscillation . 

Moore [23,24] and more recently, Hartunian and Moore 
[25] , have made preliminary studies concerned with the 
separation of the laminar boundary layer in unsteady flow. 
These studies indicate that the general case of a "semi- 
quasi-steady" oscillating flow may be divided into two 
problems, first the unsteady boundary layer perturbation 
for a fixed separation point, and second, the quasi- 
steady boundary layer over a surface which is in motion 
relative to a fixed separation point. These efforts 
failed to attain their stated goal of determining a mean- 
ingful definition of separation in unsteady flow. Hartunian 
and Moore did, however, without rigorous proof, predict the 
general effects of local accelerations and decelerations 
on the boundary layer profiles, and conjectured that 
separation might be characterized by a simultaneous vanish- 
of velocity and shear somewhere in the boundary layer. 
Perhaps the most significant feature of this work is that 
it graphically illustrates the difficulty of analytically 
arriving at a definition of separation in unsteady flow, 
as well as pointing out the need for such a definition. 

Rott [26] also comments on the problem of defining separ- 
ation in unsteady flow in a recent survey of laminar flow 
phenomena. Vidal [27] attempted to experimentally verify 
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and to extend the predictions of Hartunian and Moore. 
Unfortunately, his work was prematurely terminated, and 
although his efforts showed promise of some success, 
they must, on the whole, be considered inconclusive.. 

The phenomena of laminar boundary layer separation in 
unsteady flow remain virtually unexplored, and literally 
undefined. 
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II. LAMINAR BOUNDARY LAYER SEPARATION 



IN OSCILLATING FLOW 



A. INTRODUCTION 

To some degree, all aerodynamic flows are influenced 
by unsteadiness. Commonly encountered departures from 
steady flow cover a wide range of amplitudes and vary 
in character from the random perturbations of clear air 
turbulence to the sinusoidal oscillations encountered 
by rotary wings and turbo-machinery blades. Flow phenomena 
involving periodic unsteadiness are generally simpler to 
treat analytically and simulate experimentally than 
similar random phenomena, and any knowledge gained con- 
cerning these periodic processes will usually facilitate 
the study of their random counterparts. Consequently, 
the effects of periodic oscillations on aerodynamic flow 
phenomena are of basic interest to both researchers and 
design engineers. One of the most important of the 
phenomena affected by unsteadiness is laminar boundary 
layer separation, since its occurrence grossly alters 
the performance of flow devices. Thus, knowledge of the 
character of laminar boundary layer separation in oscillat- 
ing flow is of practical importance. 

B. CONCERNING THE DEFINITION OF SEPARATION 

It has been previously noted that little is known 
about the behavior of the laminar boundary layer in an 
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oscillating flow, adverse pressure gradient environment. 

The phenomenon of boundary layer separation in such flows 
remains virtually uninvestigated, and more particularly, 
not well defined. 

In steady flows, laminar separation is rigorously 
defined as commencing at that point on the surface where 
the first indications of fluid convection in the upstream 
direction, or "reverse flow," are observed. This reverse 
flow invariably forms eddies which decay and create the 
region of generally disturbed flow associated with 
separation. The criterion for actually determining the 
"separation point," consistent with the above definition, 
is quite logically the disappearance of shear, or equiv- 
alently, velocity gradient at the wall. 

In oscillating flows, the steady flow description of 
separation is not applicable. In such flows the character 
of the boundary layer at any spatial point varies as a 
function of time with the cyclic component of the flow. 
Hence, no single velocity profile shape can be associated 
with a given point on the surface. Moore [23] in his 
analytical investigations, predicts that instantaneous 
"zero wall shear" profiles can occur when the flow down- 
stream bears no evidence of separation. Stuart [28] calls 
attention to the fact that in oscillating flows, it is 
analytically possible to observe both instantaneous "zero 
wall shear" and "reverse flow" profiles even in the 
absence of any mean pressure gradients. This suggests that 
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regions of reverse flow occur periodically in oscillating 
boundary layers without the ensuing decay to disturbed 
flow as in the steady flow case. It is therefore obvious 
that the phenomena of laminar boundary layer separation 
must be redefined for the case of oscillating flows, and 
suitable criteria, consistent with the redefinition, must 
be determined for recognizing the occurrence of separation. 

C. OBJECTIVES 

The primary objective of the present work is to con- 
duct an experimental study to determine a rational 
definition of laminar boundary layer separation, and 
criteria for its recognition in flows wherein the free- 

/N 

stream velocity is of the form: U = U + U sin OJt. 

The method of attack consists of an investigation 
into the behavior and character of instantaneous boundary 
layer velocity profiles in such flows for as wide a 
range of freestream conditions as is practical, from which 
a rational definition and criteria should result. 

An additional objective is to correlate the occurrence 
of separation with flow parameter variation so that 
reasonable engineering predictions of separation in 
oscillating flows may be made. 
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III. EXPERIMENTAL EQUIPMENT 



A. WIND TUNNEL 

1 . General Description 

The study was conducted in the low-speed, oscillat- 
ing flow wind tunnel located in the Aeronautics Labora- 
tories of the Naval Postgraduate School. This wind tunnel 
is of open circuit design, with a 24-inch square by 223-inch 
long test section. A plan view of the tunnel is presented 
in Figure 1. The tunnel inlet is eight feet square, 
resulting in a 16:1 contraction ratio. Three high solidity 
screens located in the inlet section just upstream of the 
nozzle produce measured freestream turbulence intensities 
of from 0.261 to 0.413 per cent for the velocities encoun- 
tered in the present work. 

The wind tunnel drive consists of two Joy Axivane Fans 
in series, each of which has an internal, 100 horsepower, 
direct connected, 1750 rpm motor. The fan blades are 
internally adjustable through a pitch range of 25 to 55 
degrees, providing a wide operating base. Two sets of 
variable inlet vanes, located immediately upstream of each 
fan, are externally operated to provide control of test 
section velocity. These vanes are of multileaf design, 
and preswirl the air in the direction of fan rotation to 
reduce fan capacity. The total range of tunnel velocity 
is from 10 to 250 feet per second. 
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PLAN VIEW OF WIND TUNNEL 



2 . Rotating Shutter Valve 



Two fundamental methods of creating an oscillating 
flow environment have been employed in the past. Both 
Nickerson [7] and Hori [19] introduced oscillations by 
oscillating their models in a steady flow environment. 

This method severely restricts the range of attainable 
frequencies because of mechanical complications, and also 
introduces measurement difficulties. The other approach 
is, of course, to actually oscillate the flow over a 
stationary model. Hill [18] used a sliding shutter to 
impose oscillations on the freestream but was still 
restricted by mechanical difficulties to low frequencies. 

The most successful method of obtaining an oscillat- 
ing flow with large ranges of frequency and amplitude was 
that employed by Karlsson [29], and later by Miller [20] 
in his investigation of transition. A rotating shutter 
valve, immediately downstream of the test section, is 
employed to superimpose a periodic variation of velocity 
on the mean flow. The method used in the present investi- 
gation is virtually identical to that employed by Miller. 

The shutter valve consists of four horizontal steel shafts 
equally spaced across the test section. The shafts are 
slotted to accommodate flat blades of various widths, 
forming a set of four butterfly valves spanning the test 
section. Figure 2 is a photograph of the shutter valve. 

Each blade is driven from its immediate neighbor by means 
of a timing belt and pulley arrangement. The bottom shaft 
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FIGURE 2 

PHOTOGRAPH OF THE ROTATING SHUTTER VALVE 




is driven by a five horsepower variable speed electric 
motor, through a timing belt and pulley. An intermediate 
shaft between the motor and shutter valve permits a wide 
variety of pulley ratios. This drive arrangement provides 
a frequency range of from two cycles per second to the 
first critical frequency of 933 cycles per second. The 
electric motor presently in use, however, restricts the 
oscillation frequency to a maximum of 240 cycles per second. 
The amplitude of oscillation is controlled by blade width. 
Test section closure may be varied from 25 to 100 per cent. 
The resulting amplitude of oscillation of test section 
velocity is a function of frequency, mean velocity and 
pressure gradient. In this investigation, blades producing 
50.0, 66.7 and 82.5 per cent closure were used, resulting 
in an amplitude range of from 3 to 40 per cent of the local 
mean freestream velocity. 

3 . Test Section 

The wind tunnel test section is shown in Figure 3. 
Continuous pieces of two-inch thick aluminum, 24 inches 
wide and 223 inches long, form the upper and lower test 
section walls. Each of the side walls consists of three 
two-inch thick panels of stress-relieved Lucite. For this 
investigation, the central side wall panel on the opposite 
side of the tunnel from the control console was replaced 
with two-inch thick plywood to facilitate the mounting of 
instrumentation. The Lucite panels on the console side 
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PHOTOGRAPH OF THE TEST SECTION 



of the test section are hinged and may be raised hydraul- 
ically, providing access to the test section. The heavy 
construction of the test section is intended to minimize 
deflections induced by rapid changes in static pressure. 

Figure 4 is a typical test section velocity pro- 
file. Velocity variation is less than one per cent from 
the mean to within three inches of any wall. 

Figure 5 is an overall photographic view of the 
wind tunnel. 

B . MODEL 

The model used in this study is an elongated two- 
dimensional, symetric airfoil with a 42 inch chord and 
thickness ratio of .0357. The model is shown in Figure 6 
It was fabricated of mahogany and, with its supporting 
beams, spans the test section. From the leading edge, 
the first five inches of the model profile are a half 
ellipse terminating at the minor axis which is 1.5 inches 
long. The next sixteen inches consist of a constant 
1.5 inch thickness section. The last 21 inches consist 
of straight tapers (smoothly faired into the preceeding 
section) from 1.5 inches to a thickness of .0625 inches 
on the centerline at the trailing edge. The last twenty 
inches of the model have two three-inch wide strips of 
mild steel inlaid in the upper surface to provide a chord 
wise traversing surface for the probe carriage. These 
strips are three inches apart, centered on the mid-span 
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FIGURE 4 

TYPICAL WIND TUNNEL TEST SECTION VELOCITY PROFILE 

24 INCHES UPSTREAM OF MODEL 
UMqx = 20 ft/sec 
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FIGURE 5 

OVERALL PHOTOGRAPHIC VIEW OF THE WIND TUNNEL 
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PHOTOGRAPH OF THE MODEL 



line, and extend 1.5 inches beyond the trailing edge. 

There are seventeen static pressure taps located in the 
upper surface, distributed from leading to trailing edge 
in a line four inches from the edge of the model. 

Several features were incorporated in the model design 
to help ensure rigidity and prevent warping. The basic 
model is of laminated construction, made from .75 inch 
thick slabs running spanwise and bonded together with 
epoxy over an aluminum spar. End plates made from .0625 
inch thick aluminum T-section were incorporated, with 
the section base cut to conform to the profile shape 
and the spar imbedded in the model. The relatively 
thin trailing edge was stiffened with a four-inch wide, 
.0625 inch thick stainless steel plate inserted spanwise 
along the centerline. 

The model surface was given multiple coats of lacquer 
and hand rubbed to a mirror finish. No discontinuity 
between the inlaid steel strips and the wooden surface 
is detectable. 

For this investigation, the model was mounted hori- 
zontally on the tunnel centerline at zero angle of attack. 

Control of the pressure distribution was accomplished 
by altering the contour of the upper surface of the test 
section over the model. This was done by fastening 
shaped styrofoam blocks to the upper test section surface. 
The contour of the model surface assured the continuity 



39 



of the first and second derivatives of the chordwise 



pressure distributions thus obtained. 

C . INSTRUMENTATION 

1 . Freestream Sensors 

A conventional pitot-static tube and a fixed 
hot-wire anemometer probe were located twelve inches 
forward of the model leading edge five inches above the 
centerline and were used to determine mean freestream 
velocity. Dynamic pressure was read from a Meriam 
micro-manometer. The hot wire probe was used to record 
freestream turbulence level and to make fine adjustments 
to the mean velocity settings. 

The frequency of oscillation was obtained from 
a magnetic pickup mounted outboard of the uppermost 
shutter blade shaft, as seen in Figure 2. The output 
frequency was read on a Berkeley decade counter. 

2 . Hot-Wire Anemometers 

Two constant temperature hot-wire anemometer 
units were used in this investigation. The primary unit 
was the ten channel Security Associates Model 1000, shown 
in Figure 7. A single channel Security Associates Model 
100 hot-wire anemometer was also used as a "mobile" unit. 

It was thus possible to operate and monitor eleven inde- 
pendent hot-wire sensors simultaneously. All eleven 
basic circuits, exclusive of power supplies, potentiometers 
and switches, are identical and interchangeable solid 
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PHOTOGRAPH OF THE TEN -CHANNEL HOT WIRE ANEMOMETER 



state printed circuit plug-in modules. The circuits in 

/ 

the ten channel unit have common power supplies and 
volt-ammeter . An important feature of these units is 
that each hot-wire bridge circuit is coupled to an ana- 
logue computer circuit which inverts King's Law and pro- 
vides an output which is linear in velocity and zero for 
zero velocity. 

The "wires" used as sensing elements were 0.00015 
inch diameter tungsten filaments, .125 inches long. The 
ends of the filaments were copper plated to facilitate 
mounting on the probes (described in the next section) 
resulting in an effective sensing length of from .080 to 
.090 inches. The plating procedure and equipment is 
described in Reference [30]. 

The linearity of each of the circuits used was 
checked in steady freestream conditions using the pitot- 
static tube and freestream hot-wire probe referred to 
above. Figure 8 is a typical calibration curve. Response 
at typical boundary layer velocities was further checked 
by obtaining laminar boundary layer velocity profiles on 
the flat plate model described in Reference [31] and 
comparing the results to the Blasius profile. Figure 9 
is typical of these results. 

Since the hot-wire circuits are D.C. coupled, 
total instantaneous velocity was available at the out- 
puts. The direct current components of the outputs, 
which are proportional to the mean velocity components. 



42 



Volts 




FIGURE 8 

TYPICAL HOT WIRE ANEMOMETER CALIBRATION CURVE 

(STEADY FLOW) 



43 



V X 




u 

Uo0 



FIGURE 9 

TYPICAL HOTWIRE ANEMOMETER CALIBRATION 
VELOCITY PROFILE IN BLASIUS FLOW 
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were measured with the precision voltmeters contained in 
the two hot-wire anemometer units. The alternating current 
components, proportional to the oscillating velocity com- 
ponents, were displayed on a Tektronix 555 Dual-Beam 
oscilloscope equipped with a pair of 4-channel preamplifier 
units. Thus, of the eleven channels being monitored, 
eight of the oscillating velocity component outputs could 
be simultaneously displayed. The composition of the display 
could be easily changed by using a switching unit incor- 
porated in the output circuitry. The resulting displays 
were photographically recorded with a Tektronix oscillo- 
scope camera. 

The output signal from the shutter blade-actuated 
magnetic pickup was used as an external trigger to control 
oscilloscope sweep timing. 

Figure 10 is a schematic diagram of the instru- 
mentation circuitry. Figure 11 is a photograph of the 
control console, showing the basic instrumentation units. 

D. MULTIPLE HOT-WIRE PROBE 

1 . General Characteristics 

The prime requirement for the hot-wire probe in 
this investigation was to place a chordwise-traversable 
array of eleven hot-wires in the boundary layer. Due to 
the oscillating character of the flow environment, the 
elimination of vibration-induced relative motion between 
the hot-wires and the model surface was an important 
consideration . 
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SCHEMATIC DIAGRAM OF INSTRUMENTATION 
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PHOTOGRAPH OF THE CONTROL CONSOLE 



2 . Probe Modules 



Initial attempts at probe design centered on the 
concept of having ten hot-wire filaments, parallel to 
each other and the model surface, oriented in a single 
vertical array of a depth corresponding to expected 
boundary layer thicknesses (approximately .25 to .35 inches). 
Such an array proved to be physically impractical to con- 
struct. Even when the number of wires in the array was 
reduced to five, the mounting of the hot-wire filaments to 
the probe needles proved to be extremely difficult due to 
the proximity of the wires. In addition, these arrays 
required that the probe needles used be of such small 
diameter, and so flexible, that they would almost cer- 
tainly vibrate in response to the oscillating velocities. 

It was finally decided to divide the ten wire 
array into five two-wire subarrays extending in a spanwise 
direction. Each two-wire subarray is mounted on a separate 
plug-in probe module, as shown in Figure 12. Each module 
consists of a small block of epoxy resin, in which the 
four necessary probe needles and their associated male 
electrical contacts are cast. Jeweler's broaches were 
used for probe needles. They extend at a 45 degree angle 
relative to the module block. The tips are spaced .100 
inches apart vertically and .150 inches apart horizontally 
in a plane normal to the model surface when the module is 
mounted in the probe carriage. The hot-wire filaments were 
soldered to the probe tips, with the aid of a stereo 
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PHOTOGRAPH OF A PROBE MODULE 



microscope and specially designed holding jig, so that the 
filaments are mutually parallel. In order to eliminate 
possible displacement of the probe needles relative to 
each other, an epoxy web was formed between each two 
vertically adjacent needles. Ten of these probe modules 
were constructed and used in the investigation. 

3 . Vertical Traverse Probe 

A vertically traversable probe, controlled by an 
external micrometer through a flexible shaft arrangement, 
was also incorporated in the hot-wire array. Figure 13 
shows the traverse probe disassembled, and Figure 14, the 
micrometer control. 

4 . Probe Carriage 

The fixed array probe modules and the vertical 
traverse probe were both mounted on a chordwise-traversable 
aluminum probe carriage. Figure 15 is a photograph of 
the carriage with probes mounted. The carriage consists 
basically of an aerodynamically shaped crosspiece housing 
the probes and vertical support legs. The crosspiece is 
provided with twelve spanwise plug-in stations for the 
modules, each station having a vertical set screw which 
is used to vary the position of the hot-wires in the array 
plane. This arrangement provides a maximum fixed array 
depth of 0.40 inches normal to the model surface. The 
spanwise width of the array may be varied from 1.4 to 4.6 
inches by using various combinations of the plug-in sta- 
tions. The vertical traverse probe assembly is mounted 
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PHOTOGRAPH OF THE VERTICAL TRAVERSE PROBE DISASSEMBLED 
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FIGURE 14 

PHOTOGRAPH OF THE VERTICAL TRAVERSE MICROMETER CONTROL 
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FIGURE 15 



PHOTOGRAPH OF THE PROBE CARRIAGE 
WITH MODULES AND VERTICAL TRAVERSE PROBE 

MOUNTED 
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on the crosspiece center line such that its hot-wire always 
lies in the same plane normal to the model surface as the 
fixed array wires. 

Three teflon-coated Alnico sintered magnets are 
mounted to the outboard side of each support leg at its 
base. The magnets hold the probe carriage in positive 
contact with the steel inlays in the model, while still 
allowing the probe assembly to traverse in the chord- 
wise direction. 

The structure of the probe carriage, coupled 
with the magnetic mode of attachment to the model, virtually 
eliminates the possibility of oscillation induced relative 
motion between the probes and the model surface. 

The possibility of boundary layer profile distortion 
in the array plane due to the proximity of the probe 
carriage was investigated in a steady flow, positive 
pressure gradient environment. Figure 16 shows the results 
of this investigation. Apparent velocity distortions of 
up to five per cent were observed. Since the distortions 
were consistent across the array, and had little effect 
on the general character of the profile , they were con- 
sidered acceptable as a necessary part of the experimental 
environment. 

The probe carriage was moved fore and aft on the 
model by means of a remotely controlled traversing mech- 
anism mounted on the test section centerline aft of the 
model. The chordwise traverse range covers the aftermost 
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EFFECT OF PRESENCE OF PROBE CARRIAGE ON VELOCITY PROFILE 
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21 inches of the model. A view of the chordwise traverse 
mechanism in place is shown in Figure 17. Figure 18 shows 
the probe assembly and model mounted in the test section. 



56 




57 



PHOTOGRAPH OF THE CHORDWISE TRAVERSE MECHANISM MOUNTED IN 

THE TEST SECTION 
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FIGURE 18 

PHOTOGRAPH OF THE PROBE ASSEMBLY AND MODEL MOUNTE 

THE TEST SECTION 



IV. EXPERIMENTAL PROCEDURE 



A. CHARACTER OF THE FLOW ENVIRONMENT 

Prior to the collection of experimental data, it was 
necessary to investigate the nature of the flow field. 

The use of the spanwise distributed hot-wire array to 
observe boundary layer phenomena assumes the flow in the 
regime of interest to be two-dimensional. To determine 
if this was indeed the case, all of the array probes were 
set at the same level above the model surface. The probe 
assembly was then traversed through its full chordwise 
range for various oscillating flow conditions. This pro- 
cedure was repeated for various vextical displacements 
spanning the boundary layer and extending into the 
adjacent freestream. Comparison of output displays for 
the several wires verified that the flow was indeed two- 
dimensional throughout the experimental regime. 

Further preliminary surveys, made with the arrays of 
probes set at varying depths , revealed that the flow for 
all conditions of interest was "strictly periodic." As 
used here, "strictly periodic" refers to the fact that at 
any spatially fixed position in the flow, the wave form 
representing the oscillating component of velocity 
exhibited no aperiodic changes in its amplitude, phase, 
or shape. 

The strictly periodic quality of the flow made possible 
the employment of relatively simple, and consequently 
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accurate, methods in collecting data. While the fixed 
probe array provided data that could be used to obtain 
truly instantaneous boundary layer information, the 
relatively small number of data points available (8) , 
and the uncertainty in determining probe depth in the 
region immediately adjacent to the wall, restricts the 
detail and accuracy of the profiles obtained. However, 
since the flow was strictly periodic, the vertical 
traverse probe could be employed to provide very accurate 
profiles, that, while not in themselves instantaneous, 
are in every way equivalent to true instantaneous profiles. 
The fixed array was still of great value in qualitatively 
observing the instantaneous behavior of the boundary layer, 
as well as in supplementing and verifying observations made 
with the vertical traverse probe. 

B . INVESTIGATIVE PROCEDURE 

Preliminary investigations were carried out by 
traversing the entire adverse pressure gradient regime 
with the probe assembly. The output of the fixed probe 
array was monitored on the anemometer voltmeter and the 
oscilloscope display, providing a continuous qualitative 
picture of boundary layer behavior for a wide variety of 
flow conditions. When regions of particular interest 
were observed, accurate "instantaneous" boundary layer 
profiles were obtained by using the vertical traverse 
probe. For each point in the traverse, the steady (mean) 
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velocity component of the flow was recorded from the ane- 
mometer voltmeter, and the oscillating component was 
recorded photographically from, the oscilloscope display. 

For these photographs, the display was composed of two 
sweeps, the uppermost presenting the output of an array 
hot-wire positioned in the freestream, and the lower the 
output of the traverse probe. Boundary layer profiles 
could then be plotted by vectorially adding the mean 
velocities for each point in the traverse with the appro- 
priate oscillating component, measured from the photographs 
at the desired point in the cycle. 

C . FLOW VARIABLES 

1 . Mean Velocity 

The initial profiles acquired indicated that 
boundary layer transition effects were encountered for 
mean flow velocities of 30 feet per second and greater. 

At mean velocities less than 15 feet per second, the pro- 
files were highly susceptable to disturbances caused by 
atmospheric wind gusts at the tunnel entrance. Therefore, 
mean velocities of 15, 20 and 25 feet per second were used 
in the investigation. 

2 . Static Pressure Distribution 

Three different mean static pressure distributions 
were employed in the investigation. Figure 19 is a chord- 
wise plot of the resulting mean pressure coefficient dis- 
tributions . 
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MEAN PRESSURE COEFFICIENT DISTRIBUTIONS 
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All three distributions used were designed to pro- 
vide a significant chordwise region of favorable pressure 
gradient upstream of the beginning of the adverse pressure 
gradient region. This was done to introduce separation 
inhibiting qualities into the boundary layer, so that the 
separation process would be spatially extended, facilitat- 
ing observation of the relevent profile phenomena. The 
magnitude of the favorable pressure gradient, and its 
proximity to the adverse pressure gradient regime were 
varied in each of the three cases to provide a wide range 
of the pressure parameter. Values of the dimensionless 
pressure parameter (defined below) were graphically obtained 
from the dimensional forms of the distributions in Figure 
19. 

3 . Frequency of Oscillation 

Early in the investigation it became evident that 
for the mean velocities used, three attainable frequencies 
produced wave forms that came as close as possible to 
being sinusoidal, and were basically symmetric. These 
were 74, 142, and 226 cycles per second. 

Harmonic analyses conducted by Laferty and Hicks 
[33] on similar wave forms, produced by this equipment, 
indicate that the average contributions of second and 
third harmonics are respectively 23 and 4 per cent of the 
fundamental. Similar analyses performed by Miller [30] 
indicate that four terms of the Fourier series adequately 
represent this type waveform to within one per cent. 
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4 . Amplitude of Oscillation 



Shutter blades resulting in 87.5, 66.7, and 50.0 
per cent tunnel closure were used, producing amplitudes 
of oscillation which varied from 40 to 3 per cent mean 
frees tream velocity. 



D. DIMENSIONLESS FLOW PARAMETERS 



Commonly used forms of dimensionless velocity, fre- 
quency and amplitude were employed in this investigation, 
and are defined in the List of Symbols. The dimensionless 
pressure parameter used, is however, to the best of this 
author's knowledge, unique to this work. 

The following parameter was used to characterize the 
mean pressure gradient effects on the boundary layer: 
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The first term in the product accounts for the relative 
effects of the spatial history of the mean pressure 
gradient distribution. It is known that the presence of 
a regime of favorable pressure gradient upstream of the 
adverse pressure gradient tends to delay separation. The 
extent of this delay is increased as the magnitude of the 
favorable pressure gradient and its proximity to the 
adverse pressure gradient regime increases. With this in 
mind, the first term was formulated as the percentage of 
overall, spatially weighted, pressure gradient distribution 
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which is adverse, or separation inducing. For the special 
case of Howarth-type flow, this history term becomes equal 
to unity, due to the absence of a favorable pressure 
gradient regime. 

The complete pressure parameter accounts therefore, 
for both the magnitude of the local pressure gradient 
at separation and the influence of the preceeding pressure 
gradient distribution. 
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V. RESULTS 



A. FLOW PHENOMENA 

1 . Flow Region Classification . 

Preliminary surveys made with the multiple probe 
array revealed that the flow field produced by the adverse 
pressure gradient could be divided into three distinct 
regions. Proceeding in a streamwise direction the first 
of these regions was characterized by the regularity and 
repetitive nature of the oscillating velocity waveform 
in the boundary layer. In this region, random pertur- 
bations in the waveform were infrequent, and when observed, 
were quickly damped. When velocity fluctuations of this 
nature occurred, it was only in conjunction with similar 
perturbations in the frees tream velocity. This portion of 
the flow may be termed the "pre-wake region." The pre- 
wake region ends at a downstream point where random per- 
turbations on the oscillating velocity waveform begin to 
originate in the boundary layer, persist, and grow with 
downstream displacement. No associated velocity fluc- 
tuations are observed in the frees tream flow. This region 
will be referred to as the "region of wake formation." 
Velocity perturbations in the region of wake formation 
increase in severity until they finally begin to induce 
significant random fluctuations in the local freestream 
velocity, marking the beginning of the "region of wake 
flow. " Figure 20 presents a sequence of waveform 
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Pre-Woke 

NRe=2.745xlO; X/t=.6 1 6 




Woke Formation 



Npe =2.863x10^ X/t=.640 




Woke Formation 
NRe=2.927xl0fxA:=.655 




Wake Formation 



Nr 6 =2.991x10^ x/c=669 




Begin Woke Formation 
NRe=2.809xlC?,x/c=.626 





Woke Formotion 
NRe =2903x10; x/c=.648 




Woke Formation 

Npe =2.959x10^ x/fc=662 



Begin Wake Flow 
Nr* =3023xlC^x/fc=676 



FIGURE 20 

TYPICAL VELOCITY WAVEFORMS IN THE THREE ADVERSE 
PRESSURE GRADIENT FLOW REGIONS 
u/0 = 40 , Nf = 8.837 x 10^, Na*’.36 
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photographs along a boundary layer mean streamline, illus- 
trating the characteristics of these three distinct flow 
regions. In these photographs, as well as subsequent 
waveform sequences, the upper oscilloscope sweep depicts 
the oscillating velocity component in the local freestream. 

Also, in all waveforms presented, the voltage, or 
velocity axis is reverse-oriented, such that velocity 
increases progress in the "downward" or "minus-y" co- 
ordinate direction. This is due to voltage inversion 
inherent in the anemometer circuit. The time coordinate 
is oriented in the usual manner, i.e., increasing from 
left to right. 

2 . Flow Reversal 

Multiple probe array observations revealed that 
certain periodically occurring phenomena were consistently 
present in the boundary layer velocity waveforms. The 
most significant of these effects was characterized by an 
initial flattening of the waveform peak at the position of 
minimum local velocity. This first occurs in the lower 
half of the boundary layer. As the hot-wire is moved 
closer to the wall, the distortion intensifies, resulting 
in an apparent dip in the waveform, indicating a transient 
relative increase in velocity. This phenomenon was inter- 
preted as an indication of the existence of reverse flow 
during the affected portion of the cycle. This interpre- 
tation was based on the fact that the hot-wire anemometer 
is an absolute value sensor, hence for an instantaneous 
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reverse flow condition, the anemometer will present the 
negative, or reversed, velocity portion of the waveform as 
positive, resulting in a "folded" waveform. This deduction 
was unquestionably verified by plotting the instantaneous 
profiles for that portion of the cycle, and observing that 
the waveform dip must be considered as negative, or reversed, 
velocity to preserve the continuity of the velocity deriv- 
atives. Figure 21 presents a typical waveform sequence, 
observed during a vertical traverse of the boundary layer, 
in which the reversal phenomena is evident. 

Prior to measuring the oscillatory velocity compo- 
nents as presented in the oscilloscope photographs, the 
areas enclosed by the waveform and the calibrated reference 
line were graphically integrated. From this integration, 
corrections to the mean velocity reference position on 
the oscilloscope grid were determined and applied to the 
oscillatory velocity component measurements. The oscilla- 
tory components were measured relative to this corrected 
mean velocity reference as long as no reversal-dip was 
evident in the waveform. Since the hot-wire anemometer 
is an absolute value sensor, the occurrence of flow 
^reversal introduces an additive error in the mean velocity 
read from the anemometer voltmeter. Hence, for waveforms 
having a reversal-dip, the nodes on either side of the 
dip were recognized as establishing a reference coincident 
with zero total-instantaneous velocity, and measurements 
were then made relative to this reference, recognizing 
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y/8=l3 y/^=09 



FIGURE 21 

OSCILLATING VELOCITY COMPONENT WAVEFORMS 
RUN8|:Np|3.325xl0^Nf=3553xl^,Np=O209,N^=08 

MAJOR GRID DIVISIONS, 5MSEC X .2 VOLTS 
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that the "folded" portion of the waveform represented 
negative velocity. 

Considering anemometer accuracy and the measurement 
procedure, values of instantaneous velocity obtained in 
the above manner are estimated to have an uncertainty 
of five per cent. Due to uncertainties in the determin- 
ation of mean boundary layer thickness, the dimensionless 
normal spatial coordinate values have an estimated 
uncertainty of two per cent. 

3 . Periodic Effects 

In addition to the reversal-dip, two other period- 
ically occurring phenomena were observed in the boundary 
layer velocity waveforms. These were evident, in varying 
degrees, for most operating conditons, and are charac- 
terized as follows: 

"Amplification"-the growth, or amplification of 
the oscillatory waveform beyond its local free- 
stream magnitude as the boundary layer is 
vertically traversed; 

"Phase Shift"-a phase lead, or tendency of the 
oscillatory changes in the boundary layer to 
occur earlier in time than respective changes in 
the frees tream flow. 

Both of these effects may be observed in Figure 21. Amp- 
lification in this case is quite pronounced, and a maximum 
phase lead of 40 degrees is observed near the surface. 
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The magnitude and frequency of occurrence of these 
effects will be discussed in a subsequent section. They 
are introduced at this point because of their appearance 
in some of the results presented. 

B. INSTANTANEOUS VELOCITY PROFILES 

A total of 117 data runs were made to obtain instan- 
taneous velocity profile information. Of these, the first 
35 were made for purposes of familiarization, and to 
determine practical limits of the operating parameters. 

The remaining 82 runs yielded velocity profile sequences 
used in the investigation of separation phenomena. The 
operating parameters for these latter runs are listed 
in Table I, and representative numerical data for each run 
are presented in Appendix A. 

Figure 22 illustrates typical behavior of an instan- 
taneous velocity profile at a point in the pre-wake region, 
during a half cycle of oscillation. Figure 23 presents 
the waveform photographs from which the instantaneous 
velocity profiles were obtained. It should be noted that 
in Figure 22, as well as in sxabsequent profile sequences 
to be presented, the instantaneous boundary layer velocity 
is non-dimensionalized with respect to the local mean free- 
stream velocity, so that temporal profile variations are 
readily apparent. 

The profile sequence in Figure 22 was obtained well 
upstream of the region of wake formation. It is apparent 
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TABLE I 



OPERATING CONDITIONS FOR INSTANTANEOUS 
VELOCITY PROFILE RUNS 



RUN 




N^xlO^ 


^A 


x/c 


1 


2.745 


2.893 


.32 


.612 


2 


2.745 


5.552 


.28 


.612 


3 


2.745 


8.837 


.19 


.612 


4 


2.775 


2.893 


.25 


.619 


5 


2.775 


5.552 


.30 


.619 


6 


2.775 


8.837 


.25 


.619 


7 


2.809 


2.893 


00 

04 

• 


.626 


8 


2.809 


5.552 


.31 


.626 


9 


2.809 


8.837 


.27 


.626 


10 


2.831 


2.893 


.31 


.631 


11 


2.831 


5.552 


.31 


.631 


12 


2.831 


8.837 


.26 


.631 


13 


2.841 


2.893 


.30 


.633 


14 


2.841 


5.552 


.35 


.633 


15 


2.841 


8.837 


.30 


.633 


16 


2.155 


5.144 


.29 


.643 


17 


2.292 


5.144 


.22 


.662 


18 


2.292 


9.871 


.32 


.662 


19 


2.292 


15.709 


.22 


.662 


20 


2.767 


2.893 


00 

CN 

• 


.619 
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TABLE I (continued) 



RUN 


N„ xlO ^ 
Re 


N^xlO^ 




x/c 


21 


2.745 


2.893 


.29 


.612 


22 


2.703 


2.893 


.32 


.605 


23 


2.671 


2.893 


.32 


.598 


24 


2.639 


2.893 


.35 


.590 


25 


2.745 


5.552 


. 35 


.612 


26 


2.767 


5.552 


.37 


.619 


27 


2.799 


5.552 


.35 


.626 


28 


2.831 


5.552 


.34 


.633 


29 


2.863 


5.552 


. 33 


.640 


30 


2.895 


5.552 


. 35 


.648 


’ 31 


3.023 


8.837 


.24 


.676 


32 


2.991 


8.837 


.34 


.669 


33 


2.959 


8.837 


.36 


.662 


34 


2.927 


8.837 


.39 


.655 


35 


2.903 


8.837 


.36 


.648 


36 


2.863 


8.837 


.38 


.640 


37 


2.863 


5.552 


.36 


.640 


38 


2.852 


5.552 


.34 


.636 


39 


2.852 


5.552 


.34 


.636 


40 


2.809 


5.552 


. 36 


.626 


41 


2.075 


9.871 


.35 


.619 


42 


2.099 


9.871 


.35 


.626 


43 


2.131 


9.871 


.37 


.636 


44 


2.155 


9.871 


. 35 


.643 


45 


2.203 


9.871 


. 34 


.657 
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TABLE I (continued) 



RUN 


xlO ^ 
Re 


N^xlO^ 




x/c 


46 


2.235 


9.871 


.34 


.667 


47 


1.996 


5.144 


.32 


.595 


48 


2.020 


5.144 


.35 


.602 


49 


2.052 


5.144 


.36 


.612 


50 


2.076 


5.144 


.37 


.619 


51 


2.108 


5.144 


.34 


.629 


52 


2.148 


5.144 


.32 


.640 


53 


2.474 


15.709 


.25 


.738 


54 


2.451 


15.709 


.25 


.731 


55 


3.373 


3.553 


.28 


.610 


56 


3.341 


3.553 


.36 


.602 


57 


3.259 


1.853 


.30 


.583 


58 


3.309 


1.853 


.32 


.595 


59 


3.459 


5.655 


.30 


.619 


60 


3.373 


5.655 


.34 


.610 


61 


2.155 


5.144 


.05 


.643 


62 


2.195 


5.144 


.06 


.655 


63 


2.179 


5.144 


.06 


.650 


64 


2.219 


5.144 


.04 


.667 


65 


2.235 


5.144 


.04 


.667 


66 


2.283 


9.871 


.04 


.681 


67 


2.321 


15.709 


.03 


.690 


68 


2.826 


2.893 


.05 


.631 


69 


2.852 


5.552 


.03 


.638 


70 


2.890 


8.837 


.03 


.645 
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TABLE I (continued) 



RUN 


Nj^^xlO-5 


N^xlO^ 


^A 


x/c 


71 


3.405 


1.853 


.04 


.607 


72 


3.421 


3.553 


.04 


.612 


73 


3.453 


5.655 


.03 


.617 


74 


2.187 


5.144 


o 
1— 1 

• 


.650 


75 


2.254 


9.871 


o 
1— 1 

• 


.671 


76 


2.297 


15.709 


.08 


.686 


77 


2.801 


2.893 


.14 


.624 


78 


2.864 


5.552 


.09 


.636 


79 


2.903 


8.837 


.07 


.640 


80 


3.309 


1.853 


.09 


.593 


81 


3.325 


3.553 


00 

o 

• 


.598 


82 


3.357 


5.655 


.09 


.600 


All runs in 


Table I were 


conducted in 


Pressure 


Distribution I 
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FIGURE 22 



TYPICAL SEQUENCE OF INSTANTANEOUS 
PROFILES DURING HALF-CXLE OF 
OSCILLATION. PRE-WAKE REGION. 

WEAK TRANSIENT FLOW REVERSAL 
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y/8=IOO y/S=.84 y/^=.70 




y/^=.39 y/s=.32 y/s=26 




y/s=-24 y/s=.2i y/s=.i8 




y/^=.i3 y/fe=.07 



FIGURE 23 

OSCILLATING VELOCITY COMPONENT WAVEFORMS 

RUN47-, NRe=l996xlO^,Nf =5.I44 xI?,Na=.32 
PRE WAKE REGDN 

MAJOR GRID DIVISIONS; 5 MSEC X.2 VOLTS 
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that a reverse flow condition exists during approximately 
one-half of the full cycle of oscillation, reaching its 
maximum reversal when the instantaneous freestream velocity 
is at its minimum in the cycle. Transient flow reversal of 
this nature was found to be present to some extent in all 
adverse pressure gradients. For amplitudes of oscillation 
of less than eight per cent of the local mean freestream 
velocity, transient flow reversal first occurs at some 
point in the pre-wake region. For larger amplitudes, 
transient flow reversal was observed throughout the pre- 
wake region. This last behavior suggests that transient 
flow reversal may even occur in the zero pressure gradient 
region, as predicted by Stuart [28] and Moore [23] for 
sufficiently large amplitudes. However, being beyond the 
scope of the present work, this was not investigated. 

The profile sequence in Figure 22 illustrates a case 
of relatively weak transient flow reversal. For purposes 
of comparison. Figure 24 illustrates the nature of a strong 
transient flow reversal. Figure 25 presents the associated 

« 

waveforms. It can be seen in Figures 22 and 24 that the 
maximum and minimum velocity profiles form an envelope 
within which all profile variation is contained. 

Moving downstream, profile inflection gradually became 
more pronounced throughout the cycle. The duration of 
transient flow reversal in the cycle, however, increased 
only slightly with downstream displacement. At a point 
generally between two and four boundary layer thicknesses 
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TYRCAL SEQUENCE OF INSTANTANEOUS 
PRORLES DURING HALF-CYCLE OF 
OSCILLATION, PRE-WAKE REGION. 

STRONG TRANSIENT FLOW REVERSAL 
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y/s=55 y/s=47 y/s=.4l 




y/8=33 





y/s=.23 





y/S=l9 



y/S=06 



y/S= 15 y/ 5 =. 1 1 




y/^=.04 




FIGURE 25 

OSCILLATING VELOCITY COMPONENT WAVEFORMS 
RUN 2 ; NRe=2.745 x 10^, Nf =5.552x1^, Na=-28 
PRE WAKE REGION 

MAJOR GRID DIVISIONS; 5 MSEC X.2V0LTS 
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upstream of the region of wake formation, the reversed 
flow portion of the cycle began to expand rapidly as the 
"least reversed" profile approached a zero velocity 
gradient condition at the wall. Finally at a point within 
half a boundary layer thickness of the beginning of the 
region of wake formation, the entire cycle of profiles 
exhibited reversed flow or zero velocity gradient at the 
wall. This behavior was common to all operating conditions, 
and varied only in the magnitude of the initial cyclic 
fraction of reversed flow. Figure 26 illustrates the growth 
of the reverse flow portion of the cycle for several 
operating conditions. Figure 27 illustrates typical profile 
envelope changes with increasing transient flow reversal. 

Some mention should be made of the influence of 
amplification effects and phase leads on the instantaneous 
profile sequences. Amplification merely tends to locally 
broaden the profile envelope, without affecting the 
general profile shape, or the occurrence of transient flow 
reversal. Conversely, phase leads tend to cause a local 
narrowing of the profile envelope, in addition to that 
caused by reducing freestream amplitudes. They have little 
influence on general profile shape, even large phase leads 
change only the cyclic distribution of the transient reverse 
flow, and not its characteristic development. It may be 
generally said, that of the three prevalent waveform phe- 
nomena, only transient flow reversal significantly changes 
the character of the instantaneous velocity profile. 
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Upstream Oisplocement From RDintOf Totol Flow Reversol ^Inches 



t/j Elopsed Froction Of Freestreom Period 




RGURE 26 

EXAfy/PLES OF STREAMWISE DEVELDPEMENT OF REVERSE 
FLOW IN THE OSCILLATORY CYCLE 
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TYPICAL VARIATION OF INSTANTANEOUS PROFLE ENVELOPE 
WITH INCREASING REYNOLDS NUMBER 
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C. DEFINITION OF SEPARATION 



The importance of separation in steady flow lies in 
the fact that its occurrence, as characterized by the 
initial disappearance of shear at the wall, marks the 
beginning of a process that ultimately results in a 
streamline displacing region of wake flow. The essential 
feature of this process is the development of upstream 
fluid convection, or reverse flow, in the boundary layer, 
resulting in the production of eddies which induce wake 
formation. To be relevant, any process or phenomenon 
defined as separation in oscillating flow, must be sim- 
ilarly related to the formation of a region of wake flow. 

The findings of this investigation amply verify 
previous analytical predictions that the steady flow 
definition of separation is inapplicable in oscillating 
flows. Reverse flow, and consequently zero wall shear 
profiles, were observed to occur periodically at all 
points on the surface in the adverse pressure gradient 
regime, upstream of the region of wake formation. These 
transient flow reversals caused no observed perturbations 
either in the boundary layer or the freestream. Further, 
it can reasonably be assumed that transient flow reversal 
occurs even in the zero pressure gradient regime for 
many operating conditions. This phenomenon clearly pre- 
cludes the use of the steady flow definition, since no 
initial occurrence of a zero wall shear profile or flow 
reversal can be logically related to wake formation. 
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since time dependency is the complicating feature in 
oscillating flow, the use of a time-average or mean pro- 
file in defining separation seems logical. However, as 
has previously been discussed, mean profiles obtained with 
a hot-wire anemometer are erroneous when transient flow 
reversals are present. 

After considering all such possible alternatives, 
only one observable flow phenomenon is clearly related 
to the development of wake flow. This is the initial 
occurrence of zero velocity or reverse flow instantaneous 
profiles throughout the entire cycle of oscillation, 
or "continuous flow reversal." For all operating con- 
ditions observed, the region of wake formation invariably 
began in the immediate neighborhood of the initial point 
of continuous flow reversal. Moreover, it was found that 
this point could be rigorously and uniquely located for 
all observed operating conditions. 

The following definition is therefore proposed: 

Laminar boundary layer separation in oscillating flows 
is defined as commencing with the initial occurrence of 
zero velocity or reverse flow at some point in the velocity 
profile throughout the entire cycle of oscillation. 

This definition has the following important features: 

1. The phenomenon described has been experimentally 
observed to be directly related to the produc- 
tion of wake flow. 
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2 . 



Separation, thus defined, occurs at a unique 
point in the flow field. 

3. The definition may be extended to steady flow 
with no ambiguity. 

4. The definition is useful in analytic as well 
as experimental investigation. 

The separation point then is defined as that point where 
the first indications of time-wise continuous boundary 
layer flow reversal are observed. 

A suitable criterion for the recognition of separation 
is the observation of a velocity gradient at the wall which 
is less than or equal to zero throughout the entire cycle 
of oscillation, 

^ < 0 for y = 0 and 0 <. ± 1 

dy ~ T 

Figure 28 illustrates typical cyclic variation of the 
instantaneous velocity profile at separation. Figure 29 
presents the associated waveform photographs from which 
the profiles were derived. Figure 30 shows the change in 
character of the instantaneous velocity profile sequence 
at separation accompanied by the presence of strong flow 
reversal. Figure 31 presents the associated waveforms. 

Figures 32 through 59 show instantaneous profile 
envelopes at separation for a wide range of operating 
conditions. Figure 50 illustrates the influence of large 
phase leads on the profile envelope. These figures are 
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y/8=53 y/8=46 y/8=.40 




y/8=.36 y/8=.3l y/8=26 




y/8=23 y/8=20 y/8=.i7 




y/8=.l3 y/8=.07 

FIGURE 29 

OSCILLATING VELOCITY COMPONENT WAVEFORMS 

AT SEPARATION 

RUN 52: NRe2.l48 xlO^,Nf =5.144 xl^Np=. 1814, N^^=32 
MAJOR GRID DIVISIONS, 5 MSEC X .2 VOLTS 
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y/s=l.00 y/S=.8l y/s=.75 




y/8=67 y/s=56 y/S=50 






y/8=.22 



y/S=l4 



FIGURE 31 

OSCILLATING VELOCITY COMPONENT WAVEFORMS 

AT SEPARATION 

RUNI8:NRe=2235xlO^Nf =9.87lx I^Np=2283, 1^=32 
MAJOR GRID DIVISIONS, 5MSEC X2V0LTS 
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ENVELOPES OF INSTANTANEOUS PROFILES AT SEFARATDN 
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ENVELOPES OF INSTANTANEOUS PROFILES AT SEPARARON 
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ENVELOPES OF INSTANTANBDUS PROFILES AT SEPARATION 
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ENVELOPES OF INSTANTANEOUS PROFILES AT SEPARATION 
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ENVELOPES OF INSTANTANEOUS PROFILES AT SEPARATION 




I 



97 



ENVELOPES OF INSTANTANEOUS PROFILES AT SEPARATION 
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ENVELOPES OF WSTANTANEOUS PROFILES AT SEPARATION 
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ENVELOPES OF INSTANTANEOUS PROFILES AT SEPARATION 
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ENVELOPES OF INSTANTANEOUS PROFILES AT SERO^ATION 
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ENVELOPES OF INSTANTANEOUS PROFILES AT SEPARATION 
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ENVEUDPES OF INSTANTANEOUS PROFILES AT SEPARATON 
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ENVELOPES OF INSTANTANEOUS PROFILES AT SEPARATION 
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ENVELOPES OF NSTANTANEOUS PROFILES AT SEPARATION 
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ENVELOPE OF INSTANTANEOUS PROFILE 
AT SEFyVRATION 



representative of 27 sets of separation point data obtained 
from the original 82 profile sequences. Once the separation 
criteria were established from this information, the 
pressure distribution on the model was varied and an addi- 
tional 36 separation points located by direct reading of 
the data displays. Steady flow separation points were 
located in the usual manner for each of the combinations 
of mean velocity and pressure distribution employed in 
the investigations. These are listed in Table II. The 
change in separation point location due to the imposition 
of oscillation was then determined. All oscillating flow 
separation data are listed in Table III. 



D. CORRELATION OF SEPARATION POINT ALTERATION 

The correlation between oscillating flow separation 
point location and flow parameter variation has been 
obtained by a "least squares" fit, the details of which 
are outlined in Appendix B. 

In terms of dimensionless deviation from steady flow 
separation, the following dimensionless correlation was 
obtained: 



A 

s 



X - X 

ss s 



ss 



7 -177 

1.13 X lo' N„ ' 
Re 




N 



-.26 

f 



N 



.035 

A 



The correlation represents the data to within three per 
cent and the average scatter is less than one per cent. 

This indicates a regularity of occurrence that strongly re- 
inforces the validity of the present definition of separation. 
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TABLE II 



STEADY FLOW SEPARATION DATA 



(ft/sec) 




Np 


X /c 
ss 


Pressure Distribution I 


15 


2.427 


.3016 


. 724 


20 


3.033 


.2891 


.679 


25 


3.632 


.1083 


.650 


Pressure Distribution II 


15 


2.163 


.4485 


.645 


20 


2.746 


.3362 


.614 


25 


3.339 


.2501 


.598 


Pressure Distribution III 


15 


2.602 


.3573 


.776 


20 


3.267 


.2478 


. 731 


25 


3.885 


.0497 


.695 
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OPERATING CONDITIONS AND FLOW PARAMETERS AT SEPARATION 
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TABLE III (continued) 
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TABLE III (continued) 
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TABLE III (continued) 
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The cumulative effects of possible experimental 
inaccuracies on the correlation function were considered. 
The per cent uncertainty inherent in the correlation 
function was estimated using the method of Kline and 
McClintock [34]. The following is a brief outline of 
this procedure: 

^00 

Setting K = Np — , and rearranging the correlation 

^o 

function ; 



R=N 



-1.77 

Re 



N 



-.24 

P 



N 



-.26 

f 




_ -1.25 -1.77 -.26 1.51 -.24 -.24 _ .24 ^.035 _-.035 

= (U^) (x) (w) (V) (K) (x^) (q^) (U) (U) 



Differentiating logarithmically and dividing by R; 



^ =-1.25 ^-1.77 - 

R U X 



26 ^ + 1.51 ^ 
0) V 



-.035 f 



Taking the root square (the most probable combinatorial 
sum) ; 



9U„ 2 



T ? '^00 ^ "V ^ 

— = [(-1.25-^) +(-1.77—) +(-.26—) 

R U ^ “ 

00 



9o), ^+ 



, + (-.035^) ^ 

U 



The respective uncertainty intervals were estimated to 



be ; 



9U 



oo 




.02 





.05 



— = .002 

X 



/\ 
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9o3 

0) 



.001 



.002 



U 



.02 



9x 



Thus ; 



R 



= .034 



Hence/ there is an estimated uncertainty of 3.4 per 
cent in the correlation function's representation of the 
actual physical phenomena within the parametric limits of 
this investigation, namely; 



1.950x10^ < N„ < 3.637x10^ 
“Re 

1.852x10"^ <. < 15.709x10 

.0100 ^ Np ^ .2723 
.03 ^ ^ -44 



For this range of flow parameters the correlation 
function indicates that the oscillating flow separation 
point moves toward the steady flow value with increasing 
frequency and decreasing amplitude, as one would intui- 
tively expect. It is apparent, however, that within this 
parametric range, the effect of amplitude variation on 
separation point location is less than the estimated degree 
of uncertainty. 

Figures 59 through 62 are parametric plots of the 
correlation function. Dashed lines indicate extrapolations 
beyond the range of observed data. 
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FIGURE 59 

SERARATION POINT ALTERATION 
AS A FUNCTION OF PRESSURE PARAMETER 
AND REYNOLDS NUMBER WITH 
Nf =I.Ox»6^ 
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FIGURE 60 

SEPARATION POINT ALTERATION 
AS A FUNCTION OF PRESSURE PARAMETER 
AND REYNOLDS NUMBER WITH 
N^ = 5.0 X 10^ 
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FIGURE 61 



SEPARATION POINT ALTERATION 
AS A FUNCTION OF PRESSURE PARAMETER 
AND REYNOLDS NUMBER WITH 
Nf = 100x10^ 
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FIGURE 62 

SERARATION POINT ALTERATION 
AS A FUNCTION OF PRESSURE PARAMETER 

AND REYNOLDS NUM§ER WITH 
N =15.0x10^ 
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VI. DISCUSSION 



A. PERIODIC WAVEFORM DISTORTION 

During the course of investigation, several general 
trends in the occurrence of the two previously described 
periodic waveform effects were observed. 

1 . Amplification 

Maximum amplification at any chordwise position, 
for fixed operating conditions, generally occurred between 
the 20 per cent and 50 per cent freestream velocity 
laminae in the mean profile, similar to the findings of 
Hill [18] and Hori [19]. With all other flow parameters 
fixed, maximum amplification varied inversely with the 
amplitude of the freestream oscillation. Maximum amplifi- 
cations of 40 per cent, 100 per cent, and 200 per cent 
were attained for values of N^ in the range of .30, .10, 
and .03 respectively. In all instances of its occurrence, 
amplification was observed to increase in magnitude with 
increasing streamwise coordinate, but exhibited no identi- 
fiable trend with mean velocity or pressure gradient. 
Figure 63 presents values of maximum amplification as a 
function of freestream amplitude. Comparison with Hill's 
findings indicates that in the present work, either the 
larger adverse pressure gradient or higher frequencies, or 
both, tend to induce greater amplification in the case of 
small amplitudes. At larger amplitudes, the present 
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Moximum Amplification Of Amplitude Of Oscillotten ^ % 




FIGURE 63 

VARIATION OF MAXIMUM AMPLIFICATION WITH N^ 
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values appear to be in closer agreement with Hill's zero 
pressure gradient results. 

2 . Phase Shifts 

The phase shifts observed were invariably phase 
leads, as originally predicted by Lighthill [4], and 
observed by Nickerson [7] and Hill [18] in oscillating 
Blasius flow. For values of of the order of .30, the 
phase lead was almost non-existent, and reached a maximum 
of only 10 degrees when it did occur. Phase leads of up 
to 30 degrees were observed for values of N of approxi- 
mately .10. For values of of the order of .03, the 
maximum observed phase lead was approximately 45 degrees 
for most operating conditions. A significant exception was 

noted at these smallest amplitudes; at a Reynolds Number 

5 -5 

of .98x10 per foot and N^ of 5.144x10 , phase leads of 

up to 180 degrees were repeatedly observed, accompanied by 

amplifications of up to 200 per cent. Nothing approaching 

effects of this magnitude was observed for any other 

operating condition. Similar to amplification, phase leads 

exhibited no dependency on mean velocity or frequency, 

other than in the exceptional case mentioned above. Figure 

64 presents values of maximum phase lead as a function of 

freestream amplitude. Comparison with zero pressure 

gradient results indicate that with the exception of the 

case mentioned above, the presence of the adverse pressure 

gradient tends to suppress phase leads for freestream 

amplitudes greater than .10. 
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FIGURE 64 

VARIATION OF MAXIMUM ' PHASE LEAD WITH 



121 



B. TRANSIENT FLOW REVERSAL 



It has been generally observed, in the case of steady 
flow, that boundary layer flow reversal causes vortex 
formation, which in turn induces wake flow [32]. Observing 
the phenomena of transient flow reversal, one quite 
naturally questions the absence of disturbances in the 
boundary layer due to the reversal. One physically 
reasonable explanation is that the transient reversals 
are of such short temporal duration, that local vortices 
do not have time enough to grow to sufficient strength to 
disturb the flow. It has been observed that in the case 
of boundary layer separation on accelerating bodies , a 
finite time increment exists between initial flow reversal 
and vortex formation [32]. This line of reasoning also 
suggests that at very low frequencies , with relatively 
long periods of "transient" flow reversal, wake formation 
might conceivably begin with the initial occurrence of 
reverse flow. Such a possibility is, of course, compatible 
with the definition and theory of "quasi-steady" oscillat- 
ing boundary layer flows advanced by Lighthill [4], Lin [5] 
and others. This picture is, however, not in consonance 
with Stuart's [28] predictions of transient flow reversal 
in a zero pressure gradient, which the results of this 
investigation strongly suggest actually occur. 

Following the previous line of reasoning, this implies 
the possibility of separation occurring on a flat plate 
in low frequency, large amplitude oscillating flow, a 
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violation of the equations of motion. It seems obvious 
that much more knowledge of the interaction of time- 
dependent momentum, pressure, and viscous forces is required 
before phenomena such as transient flow reversal can be 
fully understood 
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VII. CONCLUSIONS 



From the results obtained, the following conclusions 
may be drawn: 

1. The first indications of wake formation occur at 
that streamwise point in the flow where the normal 
velocity gradient at the wall is observed, for the 
first time, to be less than or equal to zero 
throughout the entire cycle of oscillation. 

This point may then be logically defined as the 
"separation point," marking the beginning of 
boundary layer separation. 

2. The presence of an oscillating velocity component, 
superimposed on the mean flow, causes separation 
to occur upstream of the steady flow separation 
point. 

3. Increasing frequency of oscillation tends to cause 
the separation point to move downstream toward the 
steady flow location. 

4. For the range of flow parameters observed, the 
effects of amplitude variation on separation point 
location are less than the estimated experimental 
uncertainty . 

5. Transient boundary layer flow reversal occurs 
throughout the adverse pressure gradient regime 
for amplitude parameters of the order of .10 or 
greater, without inducing wake formation. 
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REPRESENTATIVE PROFILE DATA 
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20 ft/sec, 0 ) = 74 c.p.s., x = 26.0 in. from l.e., 6 = .240 
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RUN NO. 7: U =20 ft/sec, co = 74 c.p.s., x = 26.3 in. from l.e., 6 = .260 
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NO. 10: U =20 ft/sec, w = 74 c.p.s., x = 26.6 in. from l.e*/ <5 = .270 
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NO. 13: U =20 ft/sec, to = 74 c.p.s., x = 26.8 in. from l.e., 6 = .275 
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NO. 28: U = 20 ft/sec, cj = 142 c.p.s., x = 26.6 in. from l.e*/ <S = .245 











• 










• 










CN 


liO 








00 


LO 


VD 


LO 


c 


00 


lO 




CN 


O 


o 


o 


o 


•H 


o 


o 


o 


iH 


•H 


o 


o 


o 


1 — 1 


• 


• 


• 


• 




• 


• 


• 


• 




• 


• 


• 


• 






+ 


1 


LO 






+ 


1 


o 






+ 


1 










LO 










r- 












o 




lO 


CN 


LO 


o 


LO 


LO 


CN 


LO 


o 




LO 


o 


r— 1 


o 


CN 


• 


o 


1 — 1 


CN 


CN 


• 


o 


1 — 1 


fH 


iH 


# 


• 


• 


# 




• 


• 


• 


• 




• 


« 


• 


• 






+ 


1 


II 






+ 


1 


II 






+ 


1 


r- 


LO 


CN 


r- 




<5^ 


LO 


CN 


LO 






LO 


LO 


a\ 


o 


iH 


CN 


CN 




o 


iH 


CN 


CN 




o 


iH 


iH 


iH 


• 


• 


• 


• 




• 


• 


• 


# 


•>k 


• 


• 


# 


• 






+ 


1 


• 






+ 


1 


• 






+ 


1 










(U 










Q) 










00 


O 


LO 


CN 


• 


1— 1 


a\ 


o 


r- 


• 


1 — 1 


o 


a\ 


CN 


iH 


CN 


CN 


00 


iH 


1— 1 


CN 


CN 


CN 


fH 


1 — 1 


CN 


o 


CN 


• 


• 


• 


• 




* 


• 


• 


• 




• 


• 


• 


• 






+ 


1 


B 






+ 


1 


B 






+ 


1 










0 










0 










r- 


LO 


r- 


o 


u 




LO 


r- 


CD 


u 




LO 


KD 


LO 


iH 


CN 


CN 


00 


M-l 


1 — 1 


CN 


CN 


CN 


M-l 


• — 1 


CN 


i-H 


CN 


• 


• 


• 


• 




• 


• 


• 


• 




• 


• 


• 


• 






+ 


1 


• 






+ 


1 


• 






+ 


1 










c 




















CN 


o 


00 


o 


•H 


CD 


o 


o 


o 


•H 


CD 


o 


r-' 


CD 


CN 


00 


00 


00 




1 — 1 


00 


00 


00 




1 — 1 


00 


CN 


CN 


• 


• 


• 


• 


<J^ 


• 


• 


• 


• 


CN 


• 


• 


• 


• 






+ 


1 


• 






+ 


1 


• 






+ 


1 










VO 




















r- 


o 


00 




CN 


CN 


o 


o 


LO 


CN 


CN 


o 


LO 


o 


CN 




00 


00 




CN 






00 




CN 




00 


00 


• 


• 


• 


• 


II 


• 


• 


• 


• 


II 


•' • 


• 


• 


• 






+ 


1 








+ 


1 








+ 


1 










X 










X 










o 


o 


iH 


VD 




G\ 


o 


o 


1— 1 




<J\ 


o 


LO 




00 


LO 




00 




CN 


LO 


00 


00 




CN 


LO 


00 


CN 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 






+ 


1 


w 






+ 


1 


w 






+ 


1 


r- 


o 


1— 1 


o 




VD 


o 


o 


iH 


04 


LD 


o 


00 




00 


VD 






• 


00 


LD 


00 


00 


• 


00 


VD 


00 


00 


• 


• 


• 


• 


u 


• 


• 


• 


• 


u 


• 


• 


• 


• 






+ 


1 








+ 


1 








+ 


1 










CN 










CN 












o 


LO 


LO 




CN 


o 


LD 


00 




CN 


o 


r- 


o 




r- 


00 


00 


1 — 1 




r-' 


00 


00 


•H 




r-' 


00 


00 


• 


• 


• 


• 




• 


• 


• 


• 




• 


• 


• 


• 






+ 


1 


II 






+ 


1 


II 






+ 


1 


1— 1 


o 


00 




3 


C7^ 


o 


o 


o 


3 


a^ 


o 


CN 


o 


LO 


00 


00 


00 






CO 










CD 


00 


00 


• 


• 


• 


• 




• 


• 


• 


• 


•* 


• 


• 


• 


• 






+ 


1 


u 






-h 


1 


u 






+ 


1 










0) 










<u 










O 


o 


CN 


00 


to 


1— 1 


o 


LO 


LO 


CO 


iH 


o 


LO 




VD 


a\ 


00 


00 


\ 


VD 


<3^ 


00 


00 




VD 


G\ 


00 


00 


• 


• 


• 


• 




• 


• 


* 


• 


■o 


• 


• 


• 


• 






+ 


1 


m 






+ 


1 


M-l 






+ 


1 


1 — 1 


LO 


o 


CN 


o 


1— 1 


LO 


CN 


00 


o 


1 — 1 


LO 


LO 


LO 


r- 


o\ 


00 


00 


CN 


r- 


G\ 


00 


00 


CN 


r- 


G\ 


00 


00 


• 


• 


• 


• 




• 


• 


• 


• 




• 


m 


• 


• 






+ 


1 


II 






+ 


1 


II 






+ 


1 


o 


o 


CN 




. 8 


o 


o 


CN 


00 


. 8 


o 


o 


LO 


LO 


o 


o 


00 


00 


Id 


o 


o 


00 


00 


Id 


o 


o 


00 


00 


• 


• 


• 


• 




• 


• 


• 


• 




• 


• 


• 


• 


1 — i 


1 — 1 


+ 


1 


• • 


fH 


1 — 1 


+ 


1 


• • 


iH 


1 — 1 


+ 


1 




















o 


















CN 










00 














X 


c 


• 






X 


c 


• 






X 


c 






05 


-H 


o 






05 


•H 


o 






fd 


-H 






B 


e 


D 






B 


B 


D 






B 


B 


<o 


Id 


ID 


Id 


D 


‘O 


Id 


Id 


Id 


D 




Id 


Id 


Id 


\ 




\ 


\ 


D 




\ 


\ 




D 






\ 




>1 


lnJ 


C=J 


a 




>1 


I;:! 


a 


<:=! 




>1 




C=J 


<0 



134 



NO. 31: U = 20 ft/sec, w = 226 c.p.s., x = 28.1 in. from l.e., 6 = .320 
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NO. 34: U =20 ft/sec, m = 226 c.p.s., x = 27.5 in. from l.e., 6 = .280 
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RUN NO. 40: U = 20 ft/sec, o) = 142 c.p.s., x = 26.3 in. from l.e., 6 = .235 
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NO. 43; U = 15 ft/sec, u = 142 c.p.s., x = 26.7 in. from l.e., 6 = .297 
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NO. 46: U = 15 ft/sec, lo = 142 c.p.s., x = 28.0 in. from l.e., 6 = .323 
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APPENDIX B 



MATHEMATICAL OUTLINE OF PARAMETRIC CORRELATION 



The correlation between oscillating flow separation 
point location and flow parameter variation was obtained in 
the following manner: 

It was desired to determine an expression of the form: 



A 

s 



ss 



- X 



X. 



- ~ K N„^ N^^ N^^ 



ss 



Re 



N, 



( 1 ) 



A 

s 



<‘s 



K N 



Re 





(la) 



where A , N„ , N , N-, and N are experimentally 
obtained and K, a, b, c, and d are constants determined 
in such a way as to minimize the error (e) in the 
correlation . 

Simplifying the notation: 
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equation (la) becomes: 




or, in logarithmic form: 



£n X = E a. «.n y. i = 1,2,...,5 (2) 

i ^ ^ 

now considering all the separation point data: 



Z £n X . = E (E a^ £n ) i = l,2,...,5 (3) 

j ^ j i J 

j = 1,2,. ..,63 

This expression was noted to be adaptable to treatment 
by application of a modified "least squares" curve fit 
method. Accordingly, "residuals," in terms of 
observed separation point variation were defined: 



E R . = j E (E a^ £n y. . - Jin x . ) ^ 
j J ^ j i J J 



(4) 



minimizing the residual: 



k : J 



E ( E a . Jin y . . - Jin x . ) Jl y, . 
j i 1 -^13 3 n 



k = 1 , 2 , . . . , 5 



or: 



E 

j 



E a. 
1 

1 



Jin y 



13 



Jin y 



kj 



E 

j 



Jin X . 
3 



Jin y 



kj 



0 



(5) 
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in matrix notation, with: 



{A} 


= a 


5x1 


matrix of " 


a. " 
1 


[Y] 


= a 


5x63 


matrix of 


"Jin y. . 
■^13 


{X} 


= a 


6 3x1 


matrix of 


"Jin X." 
3 



equation (5) becomes; 

[Y] [Y]'^{A} - [Y] {X} = 0 (6) 

{A} = [Yy"^] ^{YX} (7) 

equation (2) becomes: 

{X} = [Y]'^{A} (8) 

and from (la) : 

{X} - {X} = {e} (9) 



Using standard matrix algebra subroutines, a computer 
solution of equation (7) was obtained, resulting in: 

^1 7 

K = e = 1.13x10 

a = -1.77 

b = -.24 

c = - . 26 

d = .035 

The desired correlation function was then: 



A 

s 



X - X 

ss s 



X 



= 1.13x10 



ss 



-1.77 .,-.24 
^Re ""P ""f 



26 .,.035 



N 



A 
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Using the solution of equation (7) , equations (8) 
and (9) were successively solved, resulting in the 
determination that the above function represents the 
data to within three per cent, with average scatter 
of only one per cent. 
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